Gravitational waves from compact binaries inspiralling along post-Newtonian accurate 

eccentric orbits: Data analysis implications 
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Compact binaries inspiralling along eccentric orbits are plausible gravitational wave (GW) sources 
for the ground-based laser interferometers. We explore the losses in the event rates incurred when 
searching for GWs from compact binaries inspiralling along post-Newtonian accurate eccentric orbits 
with certain obvious non-optimal search templates. For the present analysis, GW signals having 2.5 
post-Newtonian accurate orbital evolution are modeled following the phasing formalism, presented 
in [T. Damour, A. Gopakumar, and B. R. Iyer, Phys. Rev. D 70, 064028 (2004)]. We demonstrate 
that the search templates that model in a gauge-invariant manner GWs from compact binaries 
inspiralling under qudrupolar radiation reaction along 2PN accurate circular orbits are very efficient 
in capturing our somewhat realistic GW signals. However, three types of search templates based 
on the adiabatic, complete adiabatic and gauge-dependent complete non-adiabatic approximants, 
detailed in [P. Ajith, B. R. Iyer, C. A. K. Robinson and B. S. Sathyaprakash, Phys. Rev. D 71, 044029 
(2005)], relevant for the circular inspiral under the qudrupolar radiation reaction were found to be 
inefficient in capturing the above mentioned eccentric signal. We conclude that further investigations 
will be required to probe the ability of various types of PN accurate circular templates, employed 
to analyze the LIGO/ VIRGO data, to capture GWs from compact binaries having tiny orbital 
eccentricities. 

PACS numbers: 04.30.Db, 04.25.Nx 



I. INTRODUCTION 

Compact binaries, namely neutron star-neutron star, 
black hole-ncutron star and black hole-black hole 
binaries, inspiralling in quasi-circular orbits are the 
most plausible sources of gravitational radiation for 
the first generation ground-based laser interferometric 
gravitational- wave (GW) detectors [1-3] . GW data anal- 
ysis communities, analyzing noisy data from the operat- 
ing interferometers, require accurate and efhcient tempo- 
rally evolving GW polarizations, h^{t) and hx{t), the so 
called GW search templates. In the ongoing efforts to 
construct GW templates, inspiralling compact binaries 
are modeled as point particles moving in quasi-circular 
orbits [4]. The approximation of quasi-circularity in the 
orbital description while constructing GW templates is 
quite appropriate, because gravitational radiation reac- 
tion quickly reduces the orbital eccentricity as a compact 
binary evolves towards its last stable orbit (LSO). Em- 
ploying the dominant contributions to energy and angu- 
lar momentum losses via GW emission for an inspiralling 
compact binary, it is straightforward to deduce that when 
its semi major axis is halved, its eccentricity roughly re- 
duces by a factor of three [5]. The above argument im- 
plies that the orbital eccentricity of the Hulse- Taylor bi- 
nary pulsar when its orbital frequency reaches around 
20 Hz will be ~ 10~^. Further, it was argued that the 
templates, constructed to detect GWs from inspiralling 
compact binaries in quasi-circular orbits should be quite 
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successful in extracting GWs from (mildly) eccentric bi- 
naries [6]. Therefore, at present various GW data anal- 
ysis communities are not explicitly searching for GWs 
from compact binaries in inspiralling eccentric orbits. 

In this paper, we revisit the issue addressed in Ref. [6], 
namely the possibility of rather efficiently extracting 
GWs from compact binaries in inspiralling eccentric or- 
bits using search templates constructed for compact bi- 
naries evolving in quasi-circular orbits. In our opinion, it 
is justified to doubt the conclusions of Ref. [6] because it 
is not that difficult to note that GWs from compact bi- 
naries in non-circular orbits were not accurately modeled 
in Ref. [6], even while restricting the radiation reaction 
to the dominant order. We observed that dominant se- 
cular, but non-reactive, contributions to the orbital evo- 
lution were ignored in Ref. [6] (see our Section II for 
details). Moreover, several astrophysically motivated in- 
vestigations, carried out in the last few years, indicate 
that compact binaries in inspiralling eccentric orbits are 
plausible sources of GWs even for the ground-based GW 
detectors [7]. Therefore, it is rather important to ex- 
plore the ability of various types of circular templates to 
capture somewhat realistic modeling of GWs from com- 
pact binaries in inspiralling eccentric orbits. However, we 
would like to point out that the ability of various types 
of circular inspiral search templates, available in the LSC 
Algorithms Library (LAL) [8], to capture GW signals 
from inspiralling compact binaries having tiny orbital ec- 
centricities will require further investigations. This is due 
to the fact that in this paper, and following Ref. [6], we 
only include the dominant contributions to the reactive 
dynamics. 

This paper is organized in the following way. In Sec- 
tion II. we explicate the reasons for our current investiga- 
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tion. Section III presents a brief summary of construct- 
ing accurate gravitational waveforms for compact bina- 
ries having arbitrary mass ratio moving in inspiralling 
eccentric orbits. How we revisit the analysis, presented 
in Ref. [6], is detailed in Section IV and we also present 
our results in this section. A brief summary, our con- 
clusions and future directions are provided in Section V. 
Appendix A deals with certain computational details. 



II. RESIDUAL ORBITAL ECCENTRICITY AND 
ITS IMPLICATIONS 

The GW signals, namely and hx{t), emitted 

by inspiralling compact binaries can be accurately com- 
puted employing the post-Newtonian (PN) approxima- 
tion to general relativity. The PN approximation to the 
dynamics of inspiralling compact binaries, usually mod- 
eled to consist of point masses, provides, for example, 
the equations of motion as corrections to the Newto- 
nian one in terms of (w/c)^ ~ Gm/c^ r, where v, m and 
r are the characteristic orbital velocity, the total mass, 
and the typical orbital separation, respectively. The way 
any residual orbital eccentricity influences possible loss 
of event rate for the initial LIGO was explored, for the 
first time, in Ref. [6]. This was achieved by computing 
the drops in the signal-to-noise-ratio while searching for 
GWs from compact binaries moving in inspiralling ec- 
centric orbits with templates constructed for binaries in 
quasi-circular orbits. The search templates employed in 
Ref. [6] are essentially given by the following expressions: 



h^[M\R = -W'-^[-^\ sm20, 



(1) 



where the familiar symbols m and 77 stand for the total 
mass and the symmetric mass ratio, while C and R' de- 
note the cosine of the orbital inclination and the radial 
distance to the binary. The subscript R appearing on 
the left hand side of Eq. (1) indicates that we are using 
the so-called restricted PN waveforms as the search tem- 
plates. The temporal evolution for and w, the orbital 
phase and angular frequency, are governed by 



dt 
duj 

dt d£N I du) 



£{co) 



N 



96 



GMuj 



5/3 



OJ 



(2a) 
(2b) 



where C{uj)^ and Ej^i stand for the dominant contribu- 
tions to the GW luminosity and the Newtonian orbital 
energy and the chirp mass Ai = mrj'^^^. Further, it is 
customary to use uj = tt/qw: /gw being the GW fre- 
quency to provide the limits for uj. In the GW literature, 
Eqs. (1) and (2) provide the 'Newtonian templates' for 
the quasi-circular inspiral due to the fact that only the 
dominant contributions to C{lo) and £{uj) are required to 
construct doj/dt. However, it is interesting to note that 



Eqs. (2) provide certain 2.5PN accurate orbital phase 
evolution. This observation is due to the fact that in 
Eqs. (2), one perturbes a compact binary in an exact 
circular orbit, defined by Eq. (2a), by an expression for 
dcu/dt, given by Eq. (2b). Therefore, it is reasonable to 
state that Eqs. (1) and (2) provide a prescription to ob- 
tain 2.5PN accurate GW phasing for compact binaries 
inspiralling along exact circular orbits. We observe that 
Ref. [6] employed analytically given Fourier domain ver- 
sion of the above search templates (see Eqs. (13), (14) 
and (15) in Ref. [6]). 

Following Ref. [9], there exists two more ways of con- 
structing circular inspiral templates that incorporate the 
Newtonian reactive dynamics and we detail them below. 
The orbital phase evolution under the so-called complete 
adiabatic approximant at the Newtonian order radiation 
reaction reads 
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The crucial difference between Eq. (2) and (3) is that 
the numerator in Eq. (3b) is 2PN accurate and the use 
of the X variable: x = {Gmio/c^Y^^ . Following Ref. [9], 
we do not expand the right hand side of Eq. (3b) and 
therefore, Eqs. (3) also provide a prescription to compute 
2.5PN order orbital phase evolution. We recall that in 
the language of Ref. [9], Eq. (2) provide the adiabatic 
approximant at the Newtonian order. 

The next prescription provides the complete non- 
adiabatic approximant at the 2.5PN order and following 
[9], its orbital evolution is defined by 
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(4a) 
(4b) 

(4c) 



where y and v define the orbital separation and veloc- 
ity vectors. It is important to note that the complete 
non-adiabatic approximant is gauge-dependent and we 
employed the harmonic gauge in Eqs. (4). Further, note 
that we are using the circular limit of Eqs. (3.1) to (3.6) 
in Ref. ([9]) as we are interested in constructing only cir- 
cular templates. 
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Let us now take a closer look at the way 2.5PN accu- 
rate GW phasing for eccentric binaries, required to model 
GW signals from compact binaries in inspiralling eccen- 
tric orbits, were performed in Ref. [6]. In what follows, 
first we describe a fully PN accurate prescription to con- 
struct the above mentioned GW signals. The process of 
constructing the expected, and therefore PN accurate, 
GW signals from compact binaries inspiralling along ec- 
centric orbits, detailed in Ref. [10], is rather involved. To 
begin with, GW phasing requires an expression, similar 
to Eq. (1), for /ix- The dominant quadrupolar contribu- 
tion to /ix, denoted by /ix L,' reads 



2GmT]C 
c^R' 

2fr(j) cos 20 



+ r^(/)^ - r-" I sin2(/) 



(5) 



where f = dr/dt and = d(f>/dt. When one is inter- 
ested in doing 2.5PN accurate GW phasing for eccentric 
binaries, following Ref. [10], it can be argued that the 
dominant secular reactive evolutions for the dynamical 
variables, r, (f>, r and are governed by the following 
two differential equations: 
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where n = 2 7r/Tr and et are the so-called mean mo- 
tion and time eccentricity (and Tr being the radial or- 
bital period) and ^ is a shorthand notation for Gmn/c^. 
These orbital elements naturally arise in the PN accu- 
rate Kcplcrian type parametric solution to the conserva- 
tive PN accurate compact binary dynamics, available in 
Refs. [11-14]. However, the orbital phase of an eccentric 
binary secularly evolves partly due to the advance of pe- 
riastron and this effect appears at the first and second 
PN (conservative) orders. Therefore, in order to have 
2.5PN accurate orbital phase evolution, it is imperative 
to include secular non-radiative effects appearing at the 
first and second PN orders. In Ref. [10], this is achieved 
by having the following 2PN accurate expressions for 0, 
symbolically written as 



= A + W{l]n,et) , 
A= {l + k)n{t~ta) 
with I = n{t — to) + ci . 



c\ ■ 



(7a) 
(7b) 
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In above equations k = , A$ being the advance of pe- 
riastron in the time interval Tr- Following Ref. [10], we 
note that W{l;n,et) is a 2PN accurate function, given 
explicitly in terms of eccentric anomaly u, n and Ct and 
the mean anomaly I is related to u by the 2PN accurate 
Kepler Equation (KE). The constants to, ci and cx refer 
to some initial instant and values of I and X at t — to 



(the explicit 2PN accurate expressions for X{l;n,et) and 
W{l;n,et) in harmonic gauge is available in Eqs. (25) 
in Ref. [15]). The reactive secular changes in (f>{t) en- 
ter via n{t) and et{t) and these time dependencies are 
governed by Eqs. (6). To be consistent in a PN way, 
it is also desirable to employ 2PN accurate parameteric 
equations for r, f, (j) and in terms of n, et and u(l] n, et) 
(sec Ref. [10] for more details). Therefore to investigate 
if search templates constructed for compact binaries in 
quasi-circular orbits, given by Eqs. (1), (2), (3) and (4) 
are good enough to detect, via matched filtering, GWs 
from binaries in inspiralling eccentric orbits, it is desir- 
able to model GW signals following Ref. [10]. In other 
words, GW signals should be constructed incorporating 
2.5PN accurate orbital motion in Eq. (5). 






FIG. 1: Plots showing temporal evolution of scaled hx\Q{t), 
having t in seconds, for mi = lAMQ,m2 = WMq and ini- 
tial et = 0.1. For the upper plot, orbital motion is gov- 
erned by Eqs. (6) and (8) and we have rii = 125.7Hz and 
71/ = 1211.8Hz. For the lower plot, the orbital motion is 
2.5PN accurate and rii = 111.7Hz and n/ = 741.0Hz, obtain- 
able with the help of Eq. (9). The two sets of initial and final 
n values result from the fact that for both cases quadrupo- 
lar faw can only vary between 40Hz and (6^''^ tv G m/c'^)~^ . 
Further, notice the visible effect of k in the lower plot. 



However, this is not the way GWs from inspiralling 
eccentric orbits arc modeled in Ref. [6] and one can sum- 
marize their inaccurate construction of GW signals in the 
following way. We observe that Ref. [6] essentially em- 
ploys the same expression for /ix|q, but descriptions for 
conservative r, f, (f) and cf) were purely Newtonian and ra- 



4 



diation reaction enters via Eqs. (2) (wc recall that Rcf. [6] 
employs different orbital variables to describe the bi- 
nary dynamics and therefore required to solve numeri- 
cally three coupled differential equations). The crucial 
difference with respect to what is done in Ref. [10] is 
that the secular IPN and 2PN corrections to the orbital 
phase evolution were neglected in Ref. [6] . The conserva- 
tive orbital dynamics pursued in Rcf. [6] may be defined 
in terms of Keplerian parameterization and it reads 



"'"l^j (1-etCosw), 

. et(Gmn)i/3 . 
r =— r smu 



(8a) 



(8b) 
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(8c) 



where v — u =2 tan 



(1 — et cosm)^ ' 
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and /3 being (1 — -^/l — et )/et. In above expressions, u,v 
and et are the usual eccentric and true anomalies and 
Newtonian orbital eccentricity of the Keplerian parame- 
terization. The reactive evolution (again) is governed by 
Eqs. (6) and the explicit temporal evolution of hxlgit) 
also requires numerical solution of classical Kepler Equa- 
tion I = n{t — to) ~ u — et sinu. The above mentioned 
two approaches to obtain hxlqit) lead to quite differ- 
ent GW signals for compact binaries moving in inspi- 
ralling eccentric orbits as demonstrated in Fig. 1. For 
hxlgit), constructed following Ref. [10], one can clearly 
observe secularly changing modulations due to perias- 
tron advance. Another feature to note is that hxlqit), 
having 2.5PN accurate orbital motion, lasts longer than 
hxlq^t) having Eqs. (6) imposed on Newtonian accurate 
orbital motion. This is mainly because of using differ- 
ent Hi and Hf, initial and final values for n, for the two 
cases. This is required to make sure that the initial and 
final frequencies of the emitted quadrupolar GWs are 40 
Hz and (6'^/^ tt Gto/c^)~^ for our two cases (recall that 
advance of periastron shifts the dominant harmonic as 
explained in Ref. [16]). Further, while constructing these 
figures, we assume that first and third harmonics are ir- 
relevant for the initial LIGO. In our opinion, astrophys- 
ical GWs from compact binaries in inspiralling eccentric 
orbits are rather more accurately (and realistically) mod- 
eled in Ref. [10] compared to Ref. [6] and therefore, we are 
justified to doubt conclusions presented in Ref. [6]. How- 
ever, it is important to note that in a purely GW data 
analysis sense, while employing Eqs. (5), (6) and (8) to 
model GWs from inspiralling eccentric binaries, the con- 
clusions of Ref. [6] are indeed correct. This allowed us to 
accurately calibrate our codes with the results presented 
in Ref. [6]. 

We are now in a position to explore if the various 
circular templates, given by Eqs. (1), (2), (3) and (4), 



are going to efficient in capturing GWs from eccentric 
binaries. Because the matched filtering requires accu- 
rate secular description for GW phase evolution, this can 
roughly be done by computing and comparing (for the 
various cases under study) A/gw, the number of accu- 
mulated GW cycles between some minimum and max- 
imum GW frequencies. Making use of the definition 
■^fGW = ((/Jmax " '/>min)/7r, whcrc (/)inax and (/)min arc the 
values of the orbital phase associated with the some mini- 
mum and maximum GW frequencies, we can easily eval- 
uate the N'gw associated with an approximant in the 
time-domain. We list below the various cases involved in 
the present study, where the reactive effects are restricted 
to the dominant quadrupolar contributions. 

• Case I: The adiabatic approximant and the asso- 
ciated search templates are given by Eqs. (1) and 
(2). 

• Case II: The complete adiabatic circular inspiral 
and the associated search templates are specified 
by Eqs. (1) and (3). 

• Case HI: The gauge-dependent complete non- 
adiabatic circular inspiral and the associated search 
templates are specified by Eqs. (1) and (4). 

• Case IV: The eccentric hx{t), given by Eqs. (5), 
(8) and (6) and this is the GW signal employed in 
Ref. [6]. 

• Case V: The eccentric hx{t), given by Eqs. (5) and 
where the 2.5PN accurate orbital phase evolution 
is computed following Ref. [10]. 

• Case VI: The newly introduced circular templates 
given by Eqs. (19) and (20), 

We would like to point out that Ref. [6] investigated 
a scenario involving the cases IV and I. In this study, 
we use the case V to model the fiducial GWs from an 
eccentric binary and employ the cases I, II, HI and VI 
to construct the inspiral templates [recall that cases II, 
HI and VI do incorporate in different ways 2PN accu- 
rate conservative dynamics to construct the search tem- 
plates] . The use of the case I to model search templates in 
our comparison is based on the following two arguments. 
The first is that, strictly speaking and as explained ear- 
lier, Eqs. (2) also describe 2.5PN accurate orbital phase 
evolution of a compact binary inspiralling along exact 
circular orbits, defined by Eq. (2a). The second reason is 
that we would like to probe the performances of all the 
three types of search templates, introduced in Ref. [9] and 
having quadrupolar reactive dynamics, while treating the 
case V to model GWs from eccentric binaries. 

Let us briefly explain how we perform the Afaw com- 
putations in the time-domain for the various cases. For 

the case I, we use 40 7rHz and (6'^/^ Gm/c^) as the ini- 
tial and final values of uj for the temporal evolution de- 
fined by Eqs. (2) and similar limits also apply to the case 
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II. The difference in (j) at the above two u) values leads 
to Mow for the cases I and II. For the case III, we need 
to specify r and v values associated with the above men- 
tioned Lj limits and we use 2PN accurate expression for lj 
in terms of r, given by Eq. (3.11) in Ref. [24], to obtain 
the relevant limits for r and v = to r. For the Newto- 
nian accurate conservative and eccentric orbital motion, 
treated as the case IV, the lower and upper values of n, 
rii and n^, are identified with the above mentioned uj 
limits due to the fact that at this order n = ta. Further, 
reactive evolution is given by Eqs. (6) and the values of 
(/), obtained by using relevant parts of Eqs. (8) and the 
classical KE, evaluated at the instances when n reaches 
the values Ui and n / lead to New for the case IV. For the 
case V, the values for Ui and n j are numerically obtained 
by using the following 2PN accurate relation 
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(9) 



This is justified because of the observation in Ref. [16] 
that the dominant GW spectral component of a mildly 
eccentric compact binary, having FN accurate orbital 
motion, appears at (1 + k) n/n. Using uji ~ 40 7rHz and 

(jjf = (G^^^Gm/c^) , it is easy to obtain numerically 
Hi and Uf for the three canonical binaries considered in 
Table I. Using 2.5 FN accurate (j>{l;n,et) in harmonic 
gauge, obtainable using Eqs. (25) and (27) in Ref. [15] 
and Eqs. (6), we compute values of (f>(t) at instances 
when n reaches its above mentioned limiting values as 
0min and 4>miix- It is not very difficult to infer that we 
can choose, without any loss of generality, 0niin = at 
t = in all cases of Afcw evaluations. Finally, for the 
case VI, the limiting values for n are obtained using the 
circular limit of Eq. (9) and a procedure similar to the 
case I to compute the relevant Mgw ■ 

The resulting GW cycles, Mgw, for three typical com- 
pact binaries computed using the above mentioned six 
different prescriptions are listed in Table I. We observe 
that the differences in Mgw between cases I and IV are 
not very large. This indicates that quasi-circular tem- 
plates, given by Eqs. (1) and (2), should be highly ef- 
ficient to pick up GW signals modeled using Eqs. (5), 
(6) and (8) for small values of initial eccentricities. We 
also observe that for a given et the differences in Maw 
between the cases I and IV decrease as we increase the 
total mass. This is also consistent with the observation 
in Ref. [6] that for high mass binaries eccentricities up to 
0.2 may be tolerated while searching with quasi-circular 
templates. However, the differences in Afaw between the 
cases I, II, III and V are quite large even when initial 
et ~ 0.01. This makes it interesting to revisit the anal- 
ysis performed in Ref. [6]. However, the differences in 
Afaw between the cases V and VI indicate that the new 
type of circular inspiral templates, given by Eqs. (19) 
and (20), should be more efficient in capturing our some- 



what realistic eccentric GW signals. In the next section, 
we briefiy summarize 2.5FN accurate GW phasing, avail- 
able in Ref. [10], and present a prescription to perform it 
in a highly accurate and efficient manner. 



TABLE I: The accumulated number of GW cycles, Afaw, 
relevant for the initial LIGO, while considering the above de- 
tailed six types of orbital phase evolutions for the three types 
of canonical binaries. We recall that the radiation reaction is 
always restricted to the dominant quadrupolar order. 



TTIi/Mq : 1712 /Mq 



1.4 : 1.4 1.4 : 10 10 : 10 



Case I 



et=0 



1587.0 347.4 56.5 



Case II 



et = 



1517.7 305.8 46.6 



Case III 



et = 



1487.8 295.9 46.3 



Case IV 



et = 0.01 
et = 0.05 
et = 0.10 



1587.0 
1575.8 
1541.8 



347.4 
344.9 
337.3 



56.5 
56.1 
54.8 



Case V 



et = 0.01 
et = 0.05 
et = 0.10 



1829.5 
1817.3 
1779.2 



497.8 
494.7 
485.1 



91.9 
91.4 
89.6 



Case VI 



et = 



1830.5 497.4 91.9 



III. ACCURATE GW PHASING FOR 
ECCENTRIC BINARIES AT 2.5PN ORDER 

In what follows, we briefly summarize a prescription to 
efficiently perform 2.5PN accurate GW phasing for com- 
pact binaries of arbitrary mass ratio moving along inspi- 
ralling eccentric orbits presented in Ref. [10]. The tempo- 
ral evolution of the dynamical variables, r{t), f{t), (j){t) 
and 0(t), appearing in Eq. (5) for /ix|(5(i), is achieved 
by employing a version of the general Lagrange method 
of variation of arbitrary constants. The idea is to split 
the compact binary dynamics, defined by the relative ac- 
celeration A, into two parts. In general, A consists of 
a conservative, and in our case integrable, part and 
a reactive part A' that perturbes the conservative dy- 
namics. To perform efficient GW phasing, we employ 
a "semi-analytic" solution to the conservative dynamics 
^0- Thereafter, a solution to the PN accurate dynamics 
A — Aq -\- A! is obtained by varying the constants present 
in the generic solution to Aq. 

In the present paper, we restrict to be 2PN accurate 
and A' contains only the dominant 2.5PN contributions, 
the so called Newtonian RR terms. When the conser- 
vative binary dynamics is 2PN accurate, it is convenient 
to express the relative separation vector r (in a suitably 
defined 'center-of-mass frame') as 



(10) 
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where i = p, j = cos i q + sin i N [note that {p, q, N) is 
the same orthonormal triad employed to construct hxlq]- 
The 2PN accurate expressions for the dynamical vari- 
ables r(t), (f>{t) and their time derivatives may be ex- 
pressed parametrically as 



r{t) = r {u{l),n,et) 
f{t) = f {u{l),n,et) 



(f>{t) = X{t, 71, et) + W{u{l),n, et) 
(j)(t) = (j) {u{l),n,et) 



(11a) 
(lib) 
(11c) 
(lid) 



The explicit 2PN accurate expressions for the above 
quantities in harmonic gauge, employed in this paper, 
can be extracted from Eqs. (23), (24), (25) and (26) in 
Ref. [15] and they are explicitcly provided in appendix A. 
The basic angles I and A, symbolically given in Eq. (7) 
[see Eqs. (25) and (27) of Ref. [15] for the 2PN accurate 
explicit expressions], read 



I = n{t - to) + ci , 
X= {l + k)n{t- to) + cx 



(12a) 
(12b) 



The explicit time evolutions for r, f, (f> and (f> are provided 
by the following 2PN accurate Kepler equation, in har- 
monic gauge, which connects I and u (see Eq. (27) in 
Ref. [15]): 



I ^ u ~ et sin u 



^4/3 



(15?7 — 'rf)et sinu^J 1 — 

(13) 



8v^l — (1 — et cosm) 
+ 12(5 - 2r]){v - u){l ~ et cosu) 



In the above description for the conservative 2PN ac- 
curate dynamics, we have four constants of integration, 
namely n, et,Q and c^. When the orbital dynamics is 
fully 2.5PN accurate, i.e., ^ = + A' , then the above 
four constants of integration become time dependent. 
The differential equations governing the four constants of 
integration can be computed by demanding same func- 
tional forms for r, r and 4> even when the binary dynamics 
is fully 2.5PN accurate. Further, following Ref. [10], the 
temporal evolution in each n, et , q and cx can be treated 
to consist of a slow drift component and a rapidly oscil- 
lating part. Symbolically, the above mentioned split in 
the four variables reads 



Cait) = Ca{t) + Ca{t) 



(14) 



where a denotes one of the four constants of 2PN accu- 
rate orbital dynamics, namely, n,et,ci and cx- To our 
desired PN order, it was demonstrated in Ref. [10] that 
dci/dt = dcx/dt = and dn/dt and det/dt are given 
by Eqs. (6). For the current (and preliminary) numerical 
investigations, we neglected small amplitude fast oscilla- 
tions n^etyCx and q (however, their explicit expressions 
in harmonic gauge may be obtained from Eqs. (36) in 



Ref. [15]). To do GW phasing, we proceed as follows. 
Using 2PN accurate expressions for r, r, and (p and 2PN 
accurate Kepler equation, extractable from Eqs. (23)- 
(27) in Ref. [15], we compute r{t),r{t),(j){t) and (j){t) us- 
ing n and et to represent our PN accurate eccentric orbit. 
After that, we numerically impose reactive evolutions in 
et{t) and n{t), defined by Eqs. (6), on the 2PN accurate 
orbital dynamics and compute the associated ft-x |Q(i) us- 
ing Eq. (5). This is how we compute hxlqit) with PN 
accurate orbital evolution. 

From the above discussions, it is clear that an efficient 
implementation of GW phasing also depends on an effi- 
cient (and accurate) way of solving the 2PN accurate Ke- 
pler Equation. This is achieved by employing a slightly 
modified method of Mikkola to solve the classical Kepler 
Equation. This is detailed in the next subsection. 



A. Mikkola's solution adapted to 2PN accurate 
Kepler Equation 

There exists a plethora of analytical and numerical 
solutions associated with the celebrated classical KE, 
namely / = n{t — t^) = u — et sinu (see Ref. [17]). Ar- 
guably the most accurate and efficient numerical way to 
solve the classical KE is by Seppo Mikkola [18]. There- 
fore, we adapt Mikkola's simple and robust method to 
solve our 2PN accurate KE, given by Eq. (13). Let us 
now we briefly discuss Mikkola's procedure. Recall that 
a numerical solution to the KE usually employs Newton's 
method which requires an initial guess uq that depends 
on I and et- A number of iterations will be required to 
obtain an approximate solution that has some desired ac- 
curacy. The number of iterations to reach this accuracy 
naturally depends on uq, et and I. 

In Mikkola's method, uq is computed by introducing a 
new auxiliary variable s = sinu/3. The resulting equa- 
tion for s is Taylor expanded, forming a cubic polynomial 
in s. The resulting approximate solution is empirically 
corrected to negate the largest error occurring at / = tt 
leading to a u value which has an accuracy of not less than 
10"'^. This accuracy is greatly enhanced by employing a 
fourth order version of Newton's method (see Ref. [18] for 
further details) . This procedure provides u as a function 
of I within an accuracy of 10~^^ for all values of et and I 
for et S (0, 1). We note that the above method requires 
reduction of I into the interval — tt < ^ < tt in order to 
make s as small as possible. It is possible to employ geo- 
metrical interpretations of u and Z, to map any I in to the 
above interval, as explained in Ref. [16]. Wc note, while 
passing, that Mikkola's solution only demands the solu- 
tion of a cubic polynomial that involves no trigonometric 
functions and one-time evaluation of a couple of simple 
trigonometric functions. 

To solve the fully 2PN accurate KE, we apply 
Mikkola's method at various PN orders in a successive 
manner in the following way. First we solve the IPN ac- 
curate KE, I = n{t — to) = u — etsinw, via Mikkola's 



7 



method, and obtain u (note that due to the use of Ct 
both IPN and Newtonian accurate Kepler Equations 
have identical structure and hence solutions). This al- 
lows us to obtain 2PN corrections appearing on the right 
hand side of Eq. (13) for I in terms of (rt, e^, 1) and let us 
denote it by Z4. We obtain the solution to 2PN accurate 
KE by introducing I' = I — ~ u — et sin u and applying 
Mikkola's scheme to obtain u{l'). This allows us to get 
accurate and efficient solution to 2PN accurate KE that 
provides u in terms of I, Ct, n, m and 77. 

We are now in a position to explore data analysis impli- 
cations of hx |q(0 having 2.5PN accurate orbital motion. 



IV. NONOPTIMAL MATCHED FILTER 
SEARCHES FOR GW SIGNALS HAVING SOME 
RESIDUAL ECCENTRICITIES 

In this section, we explore the implications for ini- 
tial LIGO if hx\Q{t), given by Eq. (5), evolving under 
2.5PN accurate orbital motion, represents a potential 
GW signal. This is done by revisiting the analysis done 
in Ref. [6]. 

GW data analysts employ the technique of 'matched 
filtering' to search for GWs from compact binaries. This 
is an optimal technique provided one can construct search 
templates that accurately model expected signals from 
GW sources, especially in their phase evolution. It is 
natural that these templates depend on a number source 
parameters and therefore one needs to construct a 'bank 
of templates' that densely covers the parameter space 
of a specific (and desired) GW source. The efficiency 
of a specific template bank in detecting the associated 
GW sources is expressed by computing its Fitting Factor 
(FF), which is a measure of the 'overlap' between the 
expected GW signal and its approximate model, given 
by a template in the template bank [19]. If a template 
bank can provide FFs ^ 1, one can be confident about the 
possible detection of the associated GW signals. Further, 
it was pointed out in Rcf. [19] that the loss in event rate 
due to the deployment of inaccurate templates is cx 1 - 
FF'^. And, the typical desirable FFs for initial LIGO are 
w 0.97 [20]. 

Following Refs. [6, 19], we compute FFs in the fol- 
lowing way. Let GW signals from compact binaries in- 
spiralling along non-circular orbits be denoted by s(t) 
and let the associated quasi-circular templates be h{t; A), 
where A stands for template parameters. Further, let 
s(/) and h{f; A) denote Fourier Transforms of s{t) and 
h{t;X), computed by employing 'realft' routine of [21]. 
With these inputs and following Rcf. [19], we define the 
ambiguity function ^(A) as 



reads 



^(A) 



is\h{\)) 



(15) 



^{s\s){h{X)\h{X}) ' 
where the matched filter inner product (a|6) by definition 



ia\b) 



~a*imf) + a{f)b*if) 

Snif) 



df. (16) 



In the above expression, Snif) stands for the one sided 
noise power spectral density of the detector and we have 
for the initial LIGO 



SM) =5o[(4.49a;)-56 + 0.16a;- 
+ 0.32x2], j-^j-^^ 

= 00, f < fs, 



0.52 



(17a) 
(17b) 



where x = f/fo, with /o = 150Hz and f, = 40Hz. The 
constant Sq never enters our computations and hence 
its value is irrelevant for us. When one employs quasi- 
circular templates, given by Eqs. (1) and (2), it is possible 
to separate A into A = {to,(j)o,6), where to and (j)o de- 
note the time of arrival and the associated orbital phase 
of the binary and 6 denotes the remaining parameters. 
Following Ref. [19], one can easily maximize A over to and 
(j)Q while employing 'realft' routine to compute s(/) and 
/i(/. A). The FF is the maximum value of the ambiguity 
function, i.c, FF —"^^'^ ^(A). For circular templates, de- 
fined by Eqs. (1) and (2), the only component of 6 is the 
chirp mass A4 and therefore 



771 TTi max 



^(A). 



(18) 



However, while employing the cases II, III and IV to 
model the circular templates, it is clear that the above 
maximization has to be performed over m and 77. For 
these cases, we invoked the mini-max overlap, detailed in 
Ref. [22], and the 'amoeba' routine of Rcf. [21] to com- 
pute the associated FFs. Below we describe how wc tack- 
led the analysis performed in Ref. [6] . 



A. Brief summary of how we revisited the analysis 
of Ref. [6] 

Let us first briefly describe how GW signals were ac- 
tually constructed in Rcf. [6]. They employed Newto- 
nian accurate conservative orbital evolution and imposed 
on it the adiabatic evolution of orbital elements arising 
from Einstcinian quadrupolc formulae for the energy and 
angular momentum fluxes. The actual construction of 
GW signals from compact binaries in inspiralling eccen- 
tric orbits was achieved by employing semi-latus rectum 
p and eccentricity e to describe the Newtonian accu- 
rate (conservative) orbital motion. The effects of ra- 
diation reaction were incorporated by employing cou- 
pled differential equations for dp/dt and de/dt, com- 
puted using energy and angular momentum balances. We 
note that in Ref. [6] the construction of the GW signal, 
s{t) = hx\^{t), required solving numerically three cou- 
pled differential equations, namely, those for d(j)/dt^ dp/dt 
and de/dt (sec Eqs. (6), (8) and (9) in Rcf. [6]). Further, 
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Ref. [6] computed frequency domain versions of Newto- 
nian accurate quasi-circular templates, using the station- 
ary phase approximation [see their Eqs. (13) and (14) and 
(15)]. 

To repeat the analysis of Ref. [6] , we constructed GW 
signals by employing Eqs. (5) for hx Iqit), Eqs. (8) for the 
purely Newtonian accurate conservative orbital motion 
and Eqs. (6) to impose the effects of radiation reaction 
on the orbital motion. Therefore, our prescription leads 
to solving numerically coupled differential equations for 
n and ej along with the classical Kepler Equation. Our 
quasi-circular templates are given by Eqs. (1) and (2) 
and we employed the 'realft' routine to compute s(/) 
and h{f,X). For a given initial eccentricity, initial value 
for mean motion n = 2 tt/T is chosen following Ref. [16], 
which provided a detailed spectral analysis for GWs from 
eccentric binaries. This allowed us to choose, for Newto- 
nian binaries of all masses, rii = 125.7Hz when = 0.01, 
the initial value for et, and Ui = 83.8Hz for = 0.1. 
Note that in Ref. [6] , initial orbital frequency was always 
chosen to be 40/3 = 13.3Hz, where 40Hz being the lower 
frequency cut-off for the initial LIGO. This is mainly due 
to the assumption of Ref. [6] that GWs from eccentric bi- 
naries can be decomposed into components that oscillate 
at once, twice and thrice the orbital (radial) frequency. 
We let the final value for n be (G'^/^ Gm/c^)~^ as com- 
pact binaries that we considered here were circularized 
when n reached the above value and because n = w at 
Newtonian order. Further, in most FF computations, 
the search templates were chosen to evolve in a GW 
frequency window having 40Hz and {6^^^ n Gm/c^)~^ as 
their lower and upper limits. We noted that our results 
were somewhat insensitive to minor changes in n^. We 
are aware that Ref. [6] employed /i+ for their FF compu- 
tations. However, we choose hx mainly because we can 
factor out ~2G mC / {<?R') from our GW signals and 
templates. 

Employing the above mentioned procedure and re- 
stricting the initial values of ej < 0.1, we were able to 
reproduce most of the entries in Tables II and III of 
Ref. [6]. This gave us a lot of confidence in our routines 
to implement the FFs. In the next subsection, we detail 
the FF computations that employ Eq. (5) having 2.5PN 
accurate orbital evolution to model GWs from eccentric 
binaries, while employing the cases I, II, HI and VI as 
the search templates to probe the performances of these 
circular templates. 



B. Fitting Factor computations involving the cases 
I, II, III, V and VI 

In this section, first, we present the results from our 
FF computations that employ Eq. (5) having 2.5PN ac- 
curate orbital evolution, detailed in the Section III, to 
model, somewhat realistically, GWs from inspiralling ec- 
centric binaries. And, we employ as search templates 
the three types of circular waveforms, given by Eqs. (1) 



and (2), Eqs. (1) and (3) and Eqs. (1) and (4) [and 
these are the above detailed cases I, II and HI]. In other 
words, we are exploring if the adiabatic, complete adi- 
abatic and gauge-dependent complete non-adiabatic cir- 
cular templates, having reactive evolution restricted to 
the dominant qudrupolar order will be efficient to cap- 
ture our fiducial eccentric signal, detailed in Section HI. 
We are aware that comparing the cases V and I may be 
objectionable due to the use of Newtonian accurate quan- 
tities in the derivation of the right hand side of Eq. (2b). 
However, and as mentioned earlier, we pursued such a 
comparison due to the following two reasons. First, it is 
interesting to note, from a strict FN point of view, that 
Eqs. (2) also provide 2.5PN order GW phase evolution. 
Secondly, we want to explore whether the various types 
of circular templates, detailed in Ref. [9], can capture our 
somewhat realistic eccentric hx{t). We would like to em- 
phasize that a similar objection is not applicable while 
employing the cases II and HI to model circular tem- 
plates. The FFs, relevant for the initial-LIGO, for the 
cases I, II and HI while employing Eq. (5) having 2.5PN 
accurate orbital phase evolution to model the fiducial ec- 
centric GW signals are listed in Table II. 

Significantly lower FFs, listed in Table II, clearly indi- 
cate that GW templates, defined by Eqs. (1), (2), (3) and 
(4) are not efficient in capturing our somewhat realistic 
GW signals. For the sake of completeness, we note that 
FF computations that used GW signals having lower val- 
ues of initial et also gave us similar lower numbers [we 
note that it is rather difficult (numerically) to use values 
below 10^'^ for initial eccentricities in our GW signals 
construction]. However, low FFs are consistent with the 
arguments, based on the number of accumulated GW cy- 
cles, presented in Section II. 

Let us now explain the way we constructed the inspiral 
templates associated with the case VI. The new circular 
inspiral template family is obtained by taking the circu- 
lar limit, defined by et 0, of Eqs. (6), (26), (48d), (51), 
(52), (53) and (63) in Ref. [10]. Using these expressions 
in the limit et — > 0, we constructed a new type of GW 
templates for compact binaries inspiralling along PN ac- 
curate quasi-circular orbits. The relevant expressions are 

/ix(</',n)|fl = -4C^-2^f— sin 20, (19) 
where (\)(t) and n(i) are governed by 

| = „|:,3.- + (f-7.,)eV3|. 

In our expression for /ix (</*, 'T-)|ii: we have neglected PN 
corrections to Gm/r and (j) in the — > limit. This is 
justified as we can treat these corrections forming PN 
contributions to the amplitude of hx (this is indeed the 
way PN corrections to hx and /i+ were computed in 
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TABLE II: The FFs, relevant for the initial-LIGO, involving 
the adiabatic, complete adiabatic and gauge-dependent com- 
plete non-adiabatic approximants, namely the cases I, II and 
III, as search templates and where the reactive order is at the 
dominant quadrupolar level. The fiducial eccentric signals are 
constructed using Eq. (5) where the orbital dynamics is fully 
2.5PN accurate [the case V]. The templates that provide the 
listed FFs are characterized hy Mt, mt and rjt- We do not 
employ the third harmonic while modeling our eccentric GW 
signals. 







1.4 : 1.4 


1.4 : 10 


10 : 10 


Case V Vs Case I 


et 


= 0.01 


FF 


0.533 


0.509 


0.655 






Mt/Ms 


0.910 


0.802 


0.745 


fit 


= 0.05 


FF 


0.532 


0.507 


0.654 






Mt/Ms 


0.910 


0.802 


0.744 


et 


= 0.10 


FF 


0.522 


0.507 


0.659 






Mt/Ms 


0.912 


0.803 


0.746 


Case V Vs Case II 


et 


= 0.01 


FF 


0.430 


0.504 


0.755 






mt/rris 


0.870 


0.622 


0.614 






Vt/Vs 


0.992 


1.366 


0.987 


et 


= 0.05 


FF 


0.428 


0.500 


0.744 






mt/rris 


0.877 


0.655 


0.659 






■nt/r]s 


0.980 


1.250 


0.863 


et 


= 0.10 


FF 


0.417 


0.493 


0.732 






mt/ms 


0.900 


0.689 


0.623 






Vt/Vs 


0.939 


1.151 


0.934 


Case V Vs Case III 


et 


= 0.01 


FF 


0.424 


0.500 


0.772 






mt/ms 


0.852 


0.580 


0.601 






Vt/Vs 


1.000 


1.492 


0.997 


et 


= 0.05 


FF 


0.420 


0.494 


0.762 






mt/ms 


0.886 


0.678 


0.662 






Vt/V= 


0.936 


1.139 


0.836 


et 


= 0.10 


FF 


0.411 


0.487 


0.752 






mt/ms 


0.889 


0.737 


0.603 






Vt/Vs 


0.935 


0.984 


1.000 



Rcf. [23]). Wc detail below various propcrities of this 
new approximant. 

It is interesting to note that, in the PN approxima- 
tion and while dealing with the restricted PN waveforms. 
Eqs. (19) and (20) are simply another representation of 
Eqs. (1) and (2). In other words, both Eqs. (1) and (2) 
and Eqs. (19) and (20) can be used to describe GWs from 
quasi-circular compact binaries having 2.5PN accurate 
non-stationary orbital phase evolution. The difference is 
that the use of d(t)/dt = oj leads to a PN independent 
stationary phase evolution and while in terms of n, it 
is evidently PN dependent. Further, following [12], one 
can also demonstrate that n is indeed a gauge invari- 
ant PN accurate quantity like lo. Therefore, based on 
purely theoretical arguments, it is not possible to argue 
against the use of n to construct quasi-circular PN accu- 
rate GW search templates. In principle, one can employ 
other gauge invariant quantities like the orbital energy or 
angular momentum to construct GW search templates. 
We have used n in Eqs. (19) and (20) as 2.5PN accurate 



GW phasing for eccentric binaries, available in Ref. [10], 
employed n. And, what are the most appropriate vari- 
ables to do perform accurate GW phasing for eccentric 
binaries will be reported elsewhere. 

A rather confusing consequence of expressing d(j}/dt = 
u! = n{l + k) in terms of n for the quasi-circular inspiral 
is the following. The PN accurate conservative angular 
motion becomes (/) — (f>Q = n{l + k) {t — to)- Therefore, 
deviations from the Newtonian accurate angular motion 
can be attributed the rate of periastron advance. This is 
despite of the fact that for circular orbits, it is rather dif- 
ficult to visualize the periastron advance. However, if one 
employs, for example, the orbital energy to prescribe the 
orbit, it is no more possible to make the above statement 
as both n and k are PN accurate functions of the orbital 
energy. We note that similar arguments and expressions 
were employed in Ref. [24] during the PN accurate GW 
phasing of quasi-circular inspiral [see Eqs. (4.20)-(4.24) 
in Ref. [24] and the appendix A]. 

Finally, we would like to point out that Eqs. (20) also 
represent an adiabatic inspiral. This is because the bal- 
ance arguments, similar to the ones employed in the 
derivation of duj/dt and dx/dt, are also required to com- 
pute dn/dt. In contrast, in the complete non-adiabatic 
approximant, one does not require to invoke any balance 
arguments. Further, it is possible to provide the following 
physical explanation for the new template family. One 
may state that Eqs. (19) and (20) model GWs from a 
compact binary in a 2PN accurate circular orbit, defined 
by Eq. (20a), perturbed by the dominant order radiation 
reaction, given by Eq. (20b). We would like to recall 
that in Ref. [10], the phasing for eccentric binaries were 
performed by perturbing compact binaries in PN accu- 
rate eccentric orbits by appropriate PN accurate reactive 
dynamics for eccentric binaries. 

Therefore, we repeated our above mentioned FF com- 
putations using /i(t, A), given by Eqs. (19) and (20) and 
we again employed the 'realft' routine to compute the 
associated s{f) and h{f, A). Our results for three canon- 
ical compact binaries having three initial eccentricities, 
et = 0.01, 0.05 and 0.1 are fisted in Table. III. For mildly 
eccentric binaries, having initial et ^ 0.1, there are two 
ways of computing FFs. In the first case, we neglected 
the existence of the third harmonic which allow us to use 
Eq. (9) with lu = 40 7rHz while computing n^. Following 
Rcf. [16], we may include the third harmonic by adding 
n/2 to the right hand side of Eq. (9) and equating it to 
40 ttHz while numerically evaluating n^. The associated 
FFs are also given in Table III. 

The numbers presented in the Table III, in comparison 
with the Table. II, clearly indicate that the quasi-circular 
GW templates, given by Eqs. (19) and (20), are highly 
efficient in capturing GWs from compact binaries having 
some residual orbital eccentricity. The higher FFs for 
massive compact binaries should be due to the fact that 
they emit stronger GWs and have smaller time window 
to dephase. Finally, in Fig. 2, we plot O{0), the so-called 
overlap, obtained by maximizing the ambiguity function 
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TABLE III: The initial-LIGO FFs that uses Eqs. (19) and 
(20) to model circular GW templates and Eqs. (5) having 
2.5PN accurate orbital motion to represent fiducial eccentric 
GW signals. The symbols rrit and r]t, as expected, denote 
the total mass and the symmetric mass ratio of the desired 
template. The initial and final values of cu and n are chosen 
to have few varying between 40 Hz and {6^^^ n Gm/c^)~^ . 
For our mi = m2 — 1.4Mq binary, the orbital evolution 
is terminated when the relevant harmonic reaches lOOOHz: 
the photon shot noise limit. When the initial ej = 0.1, the 
FFs are computed by including (and neglecting) the third 
harmonic. We include the effect of the third harmonic by 
using, for example, = 78.76Hz for GWs signals having 
m = 11.4 and 77 = 0.108. 



mi/Af© : m2/M0 


1.4 : 1.4 


1.4 : 10.0 


10.0 : 10.0 


et = 0.01 FF 

mt 

m 


0.999 
2.813 
0.248 


0.999 
11.42 
0.107 


0.999 
20.16 
0.247 


et = 0.05 FF 

mt 
Vt 


0.966 
2.937 
0.233 


0.993 
11.72 
0.104 


0.992 
20.40 
0.244 


Without the third harmonic 


et = 0.10 FF 

mt 

m 


0.894 
3.133 
0.212 


0.965 
12.71 
0.094 


0.970 
20.66 
0.244 


Including the third harmonic 


et = 0.10 FF 

mt 
Vt 


0.888 
2.833 
0.246 


0.983 
11.69 
0.102 


0.984 
20.51 
0.244 



over the kinematical variables to and 0o- 

To make sure that Eqs. (19) and (20) representing our 
quasi-circular GW templates, can lead to higher FFs, we 
performed the following test. It is obvious that dn/dt, 
given in Eqs. (20), can be written in terms of 7W as 

^ = M IGMnf\2^ (21) 

Therefore, we constructed a set of quasi-circular GW 
templates that used the same m and 77 as our PN eccen- 
tric GW signal to define d(t)/dt, given in Eq. (20). How- 
ever, we used above equation for dn/dt and treated Ai as 
an independent parameter and computed FFs while max- 
imizing over M. [and this is much simpler procedure than 
maximizing over m and rf\ . The resulting FFs turned out 
to be very close to the entries listed in Table III. 

We conclude the section by providing following brief 
summary of our analysis. When one restricts radiation 
reaction to the dominant order and wants to model GWs 
from compact binaries in eccentric orbits, it is appropri- 
ate and desirable to employ Rcf. [10]. Wc have demon- 
strated that certain circular templates, given by Eqs. (19) 
and (20), representing GWs from compact binaries inspi- 
ralling under quadrupolar radiation reaction along 2PN 
accurate circular orbits, are efficient in capturing our 
fiducial GW signals from eccentric binaries. The crucial 
difference between the various circular templates, con- 
sidered as the cases I, II and III, and the case VI is that 
the Eqs. (19) and (20), representing the case VI, treat 



FIG. 2: The overlap for a binary system with mi = I.4M0, 
m2 = IO.OM0 and initial et =0.1 in the mt - r]t space. It is 
visually rather difficult to pinpoint a global maximum. How- 
ever, for a large number of {mt,Tjt) values, the overlap is more 
than 0.92. 




the conservative and reactive orbital phase evolution in 
a gauge-invariant manner with equal emphasis. 



V. CONCLUDING DISCUSSIONS 

We explored the drops in the signal-to-noise-ratios, rel- 
evant for the initial LIGO, while non-optimally search- 
ing for GWs from compact binaries inspiralling under 
quadrupolar radiation reaction along PN accurate eccen- 
tric orbits with various types of circular templates. The 
fiducial GW signals are somewhat realistically modeled 
with the help of Ref. [10], that provided /ix,+ |Q(t) from 
inspiralling eccentric binaries having 2.5PN accurate adi- 
abatic phase evolution. We demonstrated that search 
templates obtained by perturbing compact binaries in 
2PN accurate circular orbits with dominant order reac- 
tive dynamics are highly efficient in capturing our GW 
signals even from very mildly eccentric compact binaries. 
However, the search templates arising from the adiabatic, 
complete adiabatic and complete non-adiabatic approx- 
imants were found to be rather inefficient in capturing 
the fiducial GW signals having even tiny residual orbital 
eccentricities. In our view, the present analysis indicates 
that it is desirable to treat, in an equal footing, both 
the conservative and reactive contributions to the orbital 
phase evolution, while constructing the search templates 
for GWs from astrophysical compact binaries. 

We are aware that the various kinds of PN accurate 
search templates for the quasi-circular inspiral, available 
in LAL [8], employ directly or indirectly Eq. (3a) and 
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different prescriptions for the PN-accurate reactive evo- 
lution of x(t). Therefore, it will be rather important to 
extend our present analysis to investigate if various quasi- 
circular search templates, employed by the GW data an- 
alysts, can capture GW signals from compact binaries 
inspiralling along eccentric orbits and, say, having 3PN 
accurate conservative and 2PN accurate reactive dynam- 
ics. For such an analysis, the construction of PN accurate 
GW signals from inspiralling eccentric binaries should be 
influenced by Refs. [10, 15] and this is currently under ac- 
tive investigation. Further, while providing more accu- 
rate prescriptions to compute GW signals from eccentric 
binaries, it is imperative to include currently neglected 
contributions due to Cq, appearing at 2.5 and 3.5PN or- 
ders. We also feel that such post Newtonian accurate 
versions of the present analysis should be repeated for 
VIRGO, Advanced LIGO and LISA. It is also desirable 
to explore the effects of the spins on the present FF com- 
putations. This requires construction of GW signals from 
spinning compact binaries inspiralling along PN accurate 
eccentric orbits building on a preliminary investigation 
presented in Ref. [25]. 
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APPENDIX A: 2PN ACCURATE ORBITAL 
DYNAMICS FOR ECCENTRIC BINARIES 

In this appendix, we provide explicit 2PN accurate 
parametric expressions for r, r, (f) and (/) in harmonic 
gauge, required to construct hx.+ lqit) from inspiralling 
eccentric binaries having 2.5PN accurate orbital mo- 
tion. These parametric expressions are extracted from 
Eqs. (23)-(26) in Ref [15]. The radial motion, paramet- 
rically defined by r(l,n,et) and r{l,n,et), reads 



r = TN -I- riPN -I- r-2PN , 
r = rN + ?^1PN + r2PN , 



(Ala) 
(Alb) 



where 



f IPN = f N X 



(gmV'\ 

( — — I (1 - et cosw) , 



2/3 



?'2PN = m X 



6(1 — et cosu) 

^4/3 



(A2a) 

18 + 2i] - (6 - 7?7)et cos u] , (A2b) 
^^^^ _^^;;;^{-72(4-7?7)+ [72 + 30?/ + 8772 - (72 - 231?/ 4-35772 )et cos u] (1 - e^) 



72(1 - ei){l - et cosu) 



- 36(5 - 2r;)(2 4- et cos u) y^l - e^^ | , 
(G'Afn)i/3 



(1 — et cos u) 
^2/3 



et sm u , 



npN = X -^(6 - 7?/) 



?'2PN = X 



^4/3 



(A2c) 
(A2d) 

(A2e) 

468 - 15r/ + 35?/2 + (135?/ - 977^)6? + (324 4- 342r/ - 9&Tf)et cos u + (216 - 693?7 

36 



72(1 — et cosm)'^ 
1057/2) (et cosm)2 - (72 - 231?/ 4- ^WWt cosuf 



:(1 — et costz) (4 — et cos7x)(5 — 27/) 



(A2f) 



r 



The angular motion, describable using PN accurate para- 
metric expressions for and 0, is given by 
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0(A,O = A + MK(0, 

W{1) = Wn + W^lPN + W2PN 

(/> = (/>N + (/>1PN + (p2PN , 

I 



(A3a) 
(A3b) 
(A3c) 
(A3d) 



The above mentioned various PN quantities read 



k 

W2PN 



3^ 



•2/3 



^4/3 



l-e? 4(l~e?)2 
: V ~ u + Cf sin u , 
3^2/3 



[78 - 2877 + (51 - 2677)e2] 



1 



j- (w — It + 64 sin u) , 



(A4a) 
(A4b) 
(A4c) 



4/3 



32(l-e2)2(l-et cosm)3 



- p2 



</>iPN = (An X 



02PN = (p-N ^ 



(1 — et cosu)2 ' 

^2/3 



78 - 2877+(51-26?7)e2-6(5-2?7)(l-e?)3/2 (w - u)(l - e* cos u)^ (A4d) 

(A4e) 



(1 — ef){l — et cosu) 



[3 - (4 - vi)el + (1 - r])et cos u] 



(A4f) 



^4/3 



144 - 4877 - (162 + 6877 - 277^)6? + (60 + 2677 - 2077^)6^ + (I877 + 1277^)6^ 
i^ — i^iz + y (7/ — 7/ jCjj ct cos'u -|- [luo — y/7/ — 07/^ 

.2\„4 



12(l-e2)2(l-et cos7x)3 
+ [-216 + 1257; + if + (102 + 188?7 + 1677^)62 - (12 + 9777 - 77^)6^] cos it + [IO8 - 9777 - 577^ 
+ (66 - 13677 + 47/2)e2 - (48 - I777 + 17?72)e^] (e* cos uf + [-36 + 2?/ - 877^ - (6 - 7O77 - 147;2)e2] (e* cos uf 

+ 18(1 - et cosii)^(l - 2el + e* cos it) (5 - 27/)y^l-eH. (A4g) 



In the circular hmit, it is not that difficuh to conclude 
that the 2PN accurate expression for (f) becomes 

I 



i + 3e 



2/3 



39 

Y 



777 u 



c4/3 



(A5) 



The above equation, in terms of fc, is simply (p ~ 
n (1 + k). Therefore, while using 77 and in the limit 
et — > 0, deviations from the Newtonian accurate orbital 
motion may be explained in terms of k. However, this 
explanation is representation dependent as shown below. 



Using PN accurate expressions, given in Rcf. [14], one 
can easily represent in the circular limit <j) using the con- 
served orbital energy and it reads 



13 



where ( = —2E/rjm and E being the conserved center- 
of-mass orbital energy. While employing C, it is obvious 
that deviations from Newtonian accurate phase evolution 
can not be explained purely in terms of k. 

A comparison of Eqs. (A5) and (A6) leads to the fol- 
lowing observations. Using the usual definition (j) = lo = 
2ti/T and Eq. (A6), one may define a PN accurate pe- 
riod T. However, at PN orders, this is not given by 



Tr = 2 7r/n, the radial period of Keplcrian type parame- 
terization and the period relevant to construct GW sig- 
nals from inspiralling eccentric binaries. The above de- 
fined T is related to n by T = ^ (i+k) ■ Further, observe 
that there are various PN accurate ways of defining uj 
and they, at a given GW phasing order, lead to various 
quasi-circular GW templates. 
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